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CAN ONE HEAR THE SHAPE OF A DRUM? 

MARK KAC, The Rockefeller University, New York 
To George Eugene Uhlenbeck on the occasion of his sixty-fifth birthday 

“La Physique ne nous donne pas seulement 
Toccasion de resoudre des problemes . . . , elle nous 
fait presentir la solution.” H. Poincare. 

Before I explain the title and introduce the theme of the lecture I should like 
to state that my presentation will be more in the nature of a leisurely excursion 
than of an organized tour. It will not be my purpose to reach a specified des- 
tination at a scheduled time. Rather I should like to allow myself on many 
occasions the luxury of stopping and looking around. So much effort is being 
spent on streamlining mathematics and in rendering it more efficient, that a 
solitary transgression against the trend could perhaps be forgiven. 


r 



1, And now to the theme and the title. 

It has been known for well over a century that if a membrane 0, held fixed 
along its boundary T (see Fig. 1), is set in motion its displacement (in the direc- 
tion perpendicular to its original plane) 

F(x, y;t) = F(p; t) 

obeys the wave equation 

d 2 F 

= C 2 

dt 2 

where c is a certain constant depending on the physical properties of the mem- 
brane and on the tension under which the membrane is held. 

I shall choose units to make c 2 = |. 


1 
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Of special interest (both to the mathematician and to the musician) are 
solutions of the form 


Ffct) = U (p)e io>t , 

for, being harmonic in time, they represent the pure tones the membrane is capa- 
ble of producing. These special solutions are also known as normal modes. 

To find the normal modes we substitute U(p)e iwt into the wave equation 
and see that U must satisfy the equation \ V 2 £/+co 2 U — 0 with the boundary 
condition U — 0 on the boundary T of 12, corresponding to the membrane being 
held fixed along its boundary. 

The meaning of “U = 0 on T” should be made clear; for sufficiently smooth 
boundaries it simply means that U(p )—> 0 as p approaches a point of T (from the 
inside). To show that a membrane is capable of producing a discrete spectrum 
of pure tones i.e. that there is a discrete sequence of co’s coi^co 2 ^co 3 ^ • • • for 
which nontrivial solutions of 

i V 2 U + c o 2 U = 0, U = 0 on r, 


exist, was one of the great problems of 19th century mathematical physics. 
Poincare struggled with it and so did many others. 

The solution was finally achieved in the early years of our century by the 
use of the theory of integral equations. 

We now know and I shall ask you to believe me if you do not, that for regions 
12 bounded by a smooth curve T there is a sequence of numbers Xi^X 2 ^ • • • 
called eigenvalues such that to each there corresponds a function i^(p), called an 
eigenfunction, such that 

i V 2 tn + Kfn = 0 

and yp n (Jp )— >0 as p— » a point of T. 

It is customary to normalize the i^’s so that 


If, 


\f/ n (p)dp 


= 1 . 


Note that I use dp to denote the element of integration (in Cartesian co- 
ordinates, e.g., dp = dxdy). 

2. The focal point of my exposition is the following problem: 

Let 12i and 12 2 be two plane regions bounded by curves Ti and T 2 respectively, 
and consider the eigenvalue problems: 


2 V 2 C T \U — 0 in 12 1 

with 


i V 2 V + pV = 0 in 12 2 

with 


U = 0 on IT 


V = 0 on r 2 . 


Assume that for each n the eigenvalue X n for 12i is equal to the eigenvalue 
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for 12 2 . Question: Are the regions Oi and 12 2 congruent in the sense of Euclidean 
geometry? 

I first heard the problem posed this way some ten years ago from Professor 
Bochner. Much more recently, when I mentioned it to Professor Bers, he said, 
almost at once: “You mean, if you had perfect pitch could you find the shape 
of a drum.” 

You can now see that the “drum” of my title is more like a tambourine 
(which really is a membrane) and that stripped of picturesque language the 
problem is whether we can determine 12 if we know all the eigenvalues of the 
eigenvalue problem 

i \J 2 U + \U = 0 in 12, 

U = 0 on r. 

3. Before I go any further let me say that as far as I know the problem is 
still unsolved. Personally, I believe that one cannot “hear” the shape of a tam- 
bourine but I may well be wrong and I am not prepared to bet large sums 
either way. 

What I propose to do is to see how much about the shape can be inferred 
from the knowledge of all the eigenvalues, and to impress upon you the multi- 
tude of connections between our problem and various parts of mathematics and 
physics. 

It should perhaps be stated at this point that throughout the paper only 
asymptotic properties of large eigenvalues will be used. This may represent, of 
course, a serious loss of information and it may perhaps be argued that precise 
knowledge of all the eigenvalues may be sufficient to determine the shape of the 
membrane. It should also be pointed out, however, that quite recently John 
Milnor constructed two noncongruent sixteen dimensional tori whose Laplace- 
Betrami operators have exactly the same eigenvalues (see his one page note 
“Eigenvalues of the Laplace operator on certain manifolds” Proc. Nat. Acad. 
Sc., 51 (1964) 542). 

4. The first pertinent result is that one can “hear” the area of 12. This is an 
old result with a fascinating history which I shall now relate briefly. 

At the end of October 1910 the great Dutch physicist PI. A. Lorentz was 
invited to Gottingen to deliver the Wolfskehl lectures. Wolfskehl, by the way, 
endowed a prize for proving, or disproving, P'ermat’s last theorem and stipu- 
lated that in case the prize is not awarded the proceeds from the principal be 
used to invite eminent scientists to lecture at Gottingen. 

Lorentz gave five lectures under the overall title “Alte und neue Fragen der 
Physik” — Old and new problems of physics — and at the end of the fourth lecture 
he spoke as follows (in free translation from the original German): “In conclu- 
sion there is a mathematical problem which perhaps will arouse the interest of 
mathematicians who are present. It originates in the radiation theory of Jeans. 

“In an enclosure with a perfectly reflecting surface there can form standing 
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electromagnetic waves analogous to tones of an organ pipe; we shall confine our 
attention to very high overtones. Jeans asks for the energy in the frequency 
interval dv. To this end he calculates the number of overtones which lie between 
the frequencies v and v-\~dv and multiplies this number by the energy which be- 
longs to the frequency v , and which according to a theorem of statistical mechan- 
ics is the same for all frequencies. 

“It is here that there arises the mathematical problem to prove that the num- 
ber of sufficiently high overtones which lies between v and v-\-dv is independent 
of the shape of the enclosure and is simply proportional to its volume. For many 
simple shapes for which calculations can be carried out, this theorem has been 
verified in a Leiden dissertation. There is no doubt that it holds in general even 
for multiply connected regions. Similar theorems for other vibrating structures 
like membranes, air masses, etc. should also hold.” 

If one expresses this conjecture of Lorentz in terms of our membrane, it 
emerges in the form: 

iV(A) = E 1~ J — L A. 

\ n <\ 2ir 

Here 2V(X) is the number of eigenvalues less than X, |l2| the area of 12 and ^ 
means that 

N (X) | 12 1 

lim = • 

\-»<so X 2w 

There is an apocryphal report that Hilbert predicted that the theorem would 
not be proved in his life time. Well, he was wrong by many, many years. For 
less than two years later Herman Weyl, who was present at the Lorentz’ lec- 
ture and whose interest was aroused by the problem, proved the theorem in 
question, i.e. that as \-^°c 

I 

N( A) ~ L - 1 A. 

Z7T 

Weyl used in a masterly way the theory of integral equations, which his 
teacher Hilbert developed only a few years before, and his proof was a crowning 
achievement of this beautiful theory. Many subsequent developments in the 
theory of differential and integral equations (especially the work of Courant and 
his school) can be traced directly to Weyl’s memoir on the conjecture of Lorentz. 

5. Let me now consider briefly a different physical problem which too is 
closely related to the problem of the distribution of eigenvalues of the Laplacian. 

It can be taken as a basic postulate of classical statistical mechanics that if 
a system of M particles confined to a volume £2 is in equilibrium with a thermo- 
stat of temperature T the probability of finding specified particles at n, ? 2 , * * * » 

(within volume elements dr\, dr 2 , • • • , dr m) is 
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exp 


kT 


V (r h • • • , r M ) 


dr\ 


dr a 


I I 

J O a/ o 


exp 


kT 


V{r l ■ ■ ■ r M ) 


dr i • • • dr A 


where V(7 1 , • • • , 7m) is the interaction potential of the particles and k = R/N 
with R the “gas constant” and N the Avogadro number. 

For identical particles each of mass m obeying the so called Boltzmann 
statistics the corresponding assumption in quantum statistical mechanics seems 
much more complicated. One starts with the Schrodinger equation 

ti 2 2 »» ( h 

V f — V(r i, • • • , rMW = — E\j/ \ti= — > where h is the Planck constant 

2m \ 2 tt 


with the boundary condition lim • • • , 7m) =0, whenever at least one 7k 

approaches the boundary of 12. (This boundary condition has the effect of con- 
fining the particles to 12.) Let Ei^E 2 ^E^^ • • • be the eigenvalues and 
t/u, ^ 2 , • • * the corresponding normalized eigenfunctions. Then the basic postu- 
late is that the probability of finding specified particles at 7i, ? 2 , * * • , 7m (within 
dri, • • • , (Itm) is 

E e BsA TVi> • • • > r M )dn • • • dr M 

s = 1 


y e-fi s /4T 

S = 1 


There are actually no known particles obeying the Boltzmann statistics. But 
don’t let this worry you — for our purposes this regrettable fact is immaterial. 

Now, let us specialize our discussion to the case of an ideal gas which, by 
definition , means that F(?i, • • • , 7m) = 0. 

Classically, the probability of finding specified particles at n, • • ■ , r M is 
clearly 

dr i • • • drn 

where fi| is now the volume of 

Quantum mechanically the answer is not nearly so explicit. The Schrodinger 
equation for an ideal gas is 

n 2 

V 2 = - Eir 

2m 

and the equation is obviously separable. 

If I now consider the three-dimensional (rather than the 3M-dimensional) 
eigenvalue problem 
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JV¥(r) = - X^(r), 

^(r) -0 as r — > the boundary of £2, 


it is clear that the E s as well as the \[/ 8 (ri, • • ■ , n/) are easily expressible in terms 
of the X’s and corresponding ^(r)’s. 

The formula for the probability of finding specified particles at n, • • • , ym 
turns out to be 


n 


k= 1 


X) ex P 


A»» 2 ' 

mkT. 


2 


00 

X exp 

71 = 1 


X n & 2 " 
mkT. 


dr k . 


Now, as (or as T— »<x>) the quantum mechanical result should go over 

into the classical one and this immediately leads to the conjecture that as 


0 


% 2 

mkTj 


X e~ K Vn(r) 


-X n T 


/ j 

12 n= 1 


If instead of a realistic three-dimensional container 12 I consider a two- 
dimensional one, the result would still be the same 

t e~ Kr *l(?) ~ rie~ KT , r-0, 

except that now 1 12| is the area of 12 rather than the volume. 

Clearly the result is expected to hold only for r in the interior of 12. 

If we believe Weyl’s result that (in the two-dimensional case) 

I 

N(\) ~ 1 X, X — > oo , 

27T 


it follows immediately by an Abelian theorem that 


and hence that 


1 1 

— r e XnT ^ ’ r — > 0, 

12 n=l 27 TT 


n = 1 27TT 


1 

27T 



— Xr 

e 


d\. 


Setting N(X) = Xa <x^»(?), we can record the last result as 
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f* 00 1 S+ 00 

S e~ Xr dA (A) ^ — I e~ Xr d\, r — > 0. 

d o 2t d o 

Since A (A) is nondecreasing we can apply the Hardy-Littlewood-Karamata 
Tauberian theorem and conclude what everyone would be tempted to conclude, 
namely that 

A (A) = E i/v(V) ~ — ) X ^ co , 

X„<X 2t 

for every r in the interior of 12. 

Though this asymptotic formula is thus nearly “obvious” on “physical 
grounds,” it was not until 1934 that Carleman succeeded in supplying a rigorous 
proof. 

In concluding this section it may be worthwhile to say a word about the 
“strategy” of our approach. 

We are primarily interested, of course, in asymptotic properties of \ n for 
large n. This can be approached by the device of studying the Dirichlet series 

OO 

e ~ Xnt 

n = 1 

for small t. This in turn is most conveniently approached through the series 

E e Xn Vn(p) = f e~ u dA(\ ) 

n = 1 d 0 

and thus we are led to the Abelian-Tauberian interplay described above. 

6. It would seem that the physical intuition ought not only provide the 
mathematician with interesting and challenging conjectures, but also show him 
the way toward a proof and toward possible generalizations. 

The context of the theory of black body radiation or that of quantum statis- 
tical mechanics, however, is too far removed from elementary intuition and too 
full of daring and complex physical extrapolations to be of much use even in 
seeking the kind of understanding that makes a mathematician comfortable, 
let alone in pointing toward a rigorous proof. 

Fortunately, in a much more elementary context the problem of the dis- 
tribution of eigenvalues of the Laplacian becomes quite tractable. Proofs emerge 
as natural extensions of physical intuition and interesting generalizations come 
within reach. 

7. The physical context in question is that of diffusion theory , another 
branch of nineteenth century mathematical physics. 

Imagine “stuff,” initially concentrated at p( = (x 0 , yo)), diffusing through a 
plane region £2 bounded by T. Imagine furthermore that the stuff gets absorbed 
(“eaten”) at the boundary. 
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The concentration P n (p| f\ t) of matter at f{ = {x, y)) at time t obeys the 
differential equation of diffusion 


(a) 


diT 


dt 


V 2 P fi , 


the boundary condition 

(b) Pfi(p | r; t) — > 0 as r approaches a boundary point, 
and the initial condition 

(c) Pfi(p | r\ t) — > 8{r — p) as t — >0; 

here b(f — p) is the Dirac “delta function,” with “value” oc if f= p and 0 if r^p. 

The boundary condition (b) expresses the fact that the boundary is absorb- 
ing and the initial condition (c) the fact that initially all the “stuff” was con- 
centrated at p. 

I have again chosen units so as to make the diffusion constant equal to 
As is well known the concentration Po(p | r\ t) can be expressed in terms of the 
eigenvalues \ n and normalized eigenfunctions \[/ n (r) of the problem 

2 W + = 0 in 12, 

^ = 0 on r. 

In fact, P a (p\r; t ) = Y”,! e-^‘\ p a (p)rp„(r). 

Now, for small t, it appears intuitively clear that particles of the diffusing 
stuff will not have had enough time to have felt the influence of the boundary 12. 
As particles begin to diffuse they may not be aware, so to speak, of the disaster 
that awaits them when they reach the boundary. 

We may thus expect that in some approximate sense 

Pft(p | r; t) ~ P 0 (p | r; t), as t — > 0, 

where Po(p| ?; t) still satisfies the same diffusion equation 

dPo 

(a') • = i V 2 Po 


and the same initial condition 

(c') Po(p | r; t) = 8(r — p), t — > 0, 

but is otherwise unrestricted. 

Actually there is a slight additional restriction without which the solution 
is not unique (a remarkable fact discovered some years ago by D. V. Widder). 
The restriction is that P 0 ^ 0 (or more generally that P 0 be bounded from below) . 

A similar restriction for P fi is not needed since for diffusion in a bounded 
region it follows automatically. 
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An explicit formula for P 0 is, of course, well known. It is 


_ 1 

Po(p I r;t ) = — exp 
lirt 



where ||r — p|| denotes the Euclidean distance between p and r. 

I can now state a little more precisely the principle of “not feeling the bound- 
ary” explained a moment ago. 

The statement is that as t — >0 


Psi(p\r;t) = E e ut \pn(p)^n(r) 

n = 1 


2wt 


exp 



Po(p \ r; t), 


where ^ stands here for “is approximately equal to.” This is a bit vague but 
let it go at that for the moment. 

If we can trust this formula even for p = r we get 


E 


e K Yn(r) 


71=1 


1 

27 vt 


and if we display still more optimism we can integrate the above and, making 
use of the normalization condition 


obtain 



dr = 1, 


00 



71=1 


l 

2irt 


We recognize immediately the formulas discussed a while back in connection 
with the quantum-statistical-mechanical treatment of the ideal gas. If we 
apply the Hardy-Littlewood-Karamata theorem, alluded to before, we obtain 
as corollary the theorems of Carleman and Weyl. 

To do this, however, we must be allowed to interpret ^ as meaning “asymp- 
totic to.” 

8. Now, a little mathematical soul-searching. Aren’t we as far from a rigor- 
ous treatment as we were before? True, diffusion is more familiar than black 
body radiation or quantum statistics. But familiarity gives comfort, at best, 
and comfort may still be (and often is) miles away from the rigor demanded by 
mathematics. 

Let us see then what we can do about tightening the loose talk. 

First let me dispose of a few minor items which may cause you worry. 

When I write ^ = 0 on T or P(p \ r\ t )— >0 as f approaches a boundary point of 
12 there is always a question of interpretation. 
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Let me assume that T is sufficiently regular so that no ambiguity arises i.e. 
P(p | r\ t) — > 0 as r — » a boundary point of 12, 
means exactly what it says, while \p = 0 on T means 


^ — » 0 as r — > a boundary point of 12. 

Likewise, P(p\ r; /)— >5(r — p) as t— >0, has the obvious interpretation, i.e. 


lim f f PQ 
0 J J A 


r\ t) dr = 


1 


for every open set A containing p. 

Now, to more pertinent items. If the mathematical theory of diffusion corre- 
sponds in any way to physical reality we should have the inequality 


Pn(p | r;t) ^ PoCp | r;t) = 



p - r|| 2 ~ 
2 1 _ 


2wt 


For surely less stuff will be found at r at time t if there is a possibility of mat- 
ter being destroyed (on the boundary T of 12) than if there were no possibility 
of such destruction. 

Now let Q be a square with center at p totally contained in 12. Let its bound- 
ary act as an absorbing barrier and denote by P Q (p\ r\t ), the corresponding 

concentration at r at time t. 

In other words, P Q satisfies the differential equation 


(a") 


dP Q 

dt 


T v ’ p « 


and the initial condition 

(c") Pq(p \r\ t) — > S(r — p) as t — > 0. 

It also satisfies the boundary condition 

(b") P Q (p|?;/)->0 as r — » a boundary point of Q. 

Again it appears obvious that 


Pq(p \r; t) g P u (p | r; t), r E Q, 

for the diffusing stuff which reaches the boundary of Q is lost as far as P Q is con- 
cerned but need not be lost as a contribution to Pq. 

Q has been chosen so simply because Pq(p \ r\ t) is known explicitly, and, in 
particular 
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I - 4 

y 

(W 2 + % 2 ) 7T 2 

p; 0 = — 

exp < 

a 2 

m,n 

L 2a 2 J 


odd integers 


where a is the side of the square. 
The combined inequalities 


PqCp \r\ t) S P si(p ! r; t) g 


hold for all and in particular for r — p. In this case we get 


& m,n 

odd integers 


r O 2 + « 2 )tt 2 "I ~ _x t 2 r 

exp — t \ ^ 2^e 4; a {p) 


and it is a simple matter to prove that as t —> 0 we have asymptotically 
4 „ f (m 2 + n 2 ) 7r 2 "1 1 


4 T (m 2 + n 2 )Tr 2 

— 2^ ex P — — t • 

a 2 m ,n L 2a 2 


odd integers 


Thus asymptotically for t-^0 ^” =1 e~ Xnt ^\{p) / litt and Carleman’s theorem 

follows. 

It is only a little harder to prove Weyl’s theorem. 

If one integrates over Q the inequality 


m ex p - 

m,n — 


(w 2 + « 2 )ir 2 1 “ —\ n t 2 ,~. 

— 1 ^ Up) 


one obtains 


4 E exp 


( m 2 + n 2 )Tr 2 
2a? 1 


= E e X " f f 'f'l(p) dp- 

n= 1 J J Q 


We now cover 12 with a net of squares of side a , as shown in Fig. 2, and keep 
only those contained in 12. Let N(a) be the number of these squares and let 
12(a) be the union of all these squares. We have 

oo oo n> oo /*> n> 

E e- x ” ! = E VCp) dp ^ y e-^ I **(?) dp 

n—1 n = 1 J J SI n = 1 J ^Q(a) 

„ r (m 2 + n 2 ) 7T 2 
^ 4 N(a) y exp — 1 

m.n L 2a _ 
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and, integrating the inequality P Q (p\p; t)^l/2irt over 12 we get i e~* nt 
^ 1 12| /2irt. 



Noting that N(a)a 2 = 1 12(a) | we record the fruits of our latest labor in the 
form of the inequality 


G(a) — 2^ ex P 


n 2 

W' m,n 


(m 2 + w 2 )tt 2 i 


2a 2 


= e~^ nt ^ 


— 1 n=l 


I 

27r/ 


From the fact (already noted above) that 

(in 2 + n 2 )7r 2 


lim 2 irt 
t-* o 


4 v r 

— L exp - 

& m,n L 


2 a 2 


= 1 


odd 


we conclude easily that 

00 CO 

| 12(a) | rg lim inf 2irt e~ Xnt rg lim sup 2 irt ^ e -Xn * =* | 0 | ; 

*0 n — 1 t — »0 7i=l 

and since, by choosing a sufficiently small, we can make 1 12(a) | arbitrarily close 
to 1 12 1 , we must have lim^ 0 2 irt e~* nt = 1 12 1 or, in other words, 


p — ^ nt , 


12 


t-> 0. 


7i= i 2irt 

9. Are we now through with rigor? Not quite. For while the inequalities 


ex P 


r-p 


PaG | r; t) g 


2< J 


2irt 
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Ai(p| r\ t) ^ PqCp\ r; t), r E Q, 

are utterly obvious on intuitive grounds they must be proved. Let me indicate 
a way of doing it which is probably by far not the simplest. I am choosing it to 
exhibit yet another physical context. 

It has been known since the early days of this century, through the work of 
Einstein and Smoluchowski, that diffusion is but a macroscopic manifestation 
of microscopic Brownian motion. 

Under suitable physical assumptions Pq(p\ r; t) can be interpreted as the 
probability density of finding a free Brownian particle at r at time t if it started 
on its erratic journey at / = 0 from p and if it gets absorbed when it comes to the 
boundary of 12. 

If a large number N of independent free Brownian particles are started from 
p then 

N fL P(p | r ; t) dr 


is the average number of these particles which are found in A at time t. Since 
the statistical percentage error is of the order 1 /y/N continuous diffusion theory 
is an excellent approximation when N is large. 

A significant deepening of this point of view was achieved in the early 
twenties by Norbert Wiener. Instead of viewing the problem as a problem in 
statistics of particles he viewed it as a problem in statistics of paths . Without 
entering into details let me review briefly what is involved here. 

Consider the set of all continuous curves r(r), O^r < co , starting from some 
arbitrarily chosen origin 0. Let 12i, 12 2 , • • • , Q n be open sets and t\ <t 2 < •••<£„ 
ordered instants of time. The Einstein-Smoluchowski theory required that 
(with suitable units) 


Prob. {p + r(ti) G 12i, p + r{t^) G 12 2 , • • • , p + r(t n ) G 12n} 

= 1 • • • I A (p I ri;h)P 0 (ri \ r 2 ; h — h) ■ ■ ■ Po(r n - 1 1 r„\t n — t n -i) dr i • • • dr n 

do, J O 


where, as before, 


-i - 1 

A(p | r; t) = — exp 
Zirt 



Wiener has shown that it is possible to construct a completely additive mea- 
sure on the space of all continuous curves f(r) emanating from the origin such 
that the set of curves p + r(r) which at times £i<£ 2 < • • • <t n find themselves 
in open sets S2i, £2 2 , ■ • • , respectively, has measure given by the Einstein- 
Smoluchowski formula above. 

The set of curves such that p + r(r)G^, and p + ?(£)G^4 (A — an 

open set) turns out to be measurable and it can be shown, if 12 has sufficiently 
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smooth boundaries, that this measure is equal to 

I -Pn(p | r; t) dr. 

J A 

This is not a trivial statement and it should come as no surprise that it 
trivially implies the inequalities we needed a while back to make precise the 
principle of not feeling the boundary. 

In fact, as the reader no doubt sees, the inequalities in question are simply a 
consequence of the fact that if sets Cfc, (B, (3 are such that 

Ct C (B C <3 

then meas. Cf^meas. (B^meas. (3. 

One final remark before we go on. The set of curves for which 

P + r(r) G 0 ^ r ^ t and p + r(/) G A 

is measurable even if the boundary of 12 is quite wild. The measure can still be 
written as JaPq(p | f\ t)dr and it can be shown that in the interior of 12, P n (p | ?; t) 
satisfies the diffusion equation dPn/dt = l V 2 Po as well as the initial condition 

lim I Pfl(p; ?; /) dr = 1, 

*->o 

for all open sets A such that pG^4. 

It is, however, no longer clear how to interpret the boundary condition that 
Po(p | r; t) — » 0 when r — > T. 

This difficulty forces the classical theory of diffusion to consider reasonably 
smooth boundaries. The probabilistic interpretation of P n (p| r; t) provides a 
natural definition of a generalized solution of the boundary value problem under 
consideration. 

10. We are now sure that we can hear the area of a drum and it may seem 
that we spent a lot of effort to achieve so little. 

Let me now show you that the approach we used can be extended to yield 
more, but to avoid certain purely geometrical complications I shall restrict my- 
self to convex drums. 

We have achieved our first success by introducing the principle of not feeling 
the boundary. But if p is close to the boundary Y of 12 then the diffusing particles 
starting from p will, to some extent, begin to be influenced by T. 

Let q be the point on T closest to p and let Z(p) be the straight line perpendicu- 
lar to the line joining p and q. (See Fig. 3.) Then a diffusing particle starting from 
p will see for a short time the boundary T as the straight line /(p). 

One may say, using again somewhat picturesque language, that, for small t, 
the particle has not had time to feel the curvature of the boundary. 
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If this principle is valid I should be allowed to approximate (for small t) 

Pn(p\r;t) by PiC P )(p \r\t), 

where Pic p )(p \?\ 0 satisfies again the diffusion equation 

dP 1 

= — V 2 P 

dt 2 


with the initial condition P— >5(p — r) as >0, but with the boundary condition 
Pi(p)(p | r\ t) — > 0 as r approaches a point on /(p). 




Carrying this optimism as far as possible we would expect that to a good 
approximation 


/ 

J f] 


L 


Pq(p \ p',t) dp~ I PiCp)(p | p; 0 dp 


It is well known that 


PlCp)(P I P'yt) = 


1 - e~ ui l‘ 

2wt 


where 5= q — p = minimal distance from p to T. Thus (hopefully!) 


f f Pa(p | p;t) dp = E < 

J J n n—1 


I 12 | 1 

2i rt 2irt 


L 


-28*/t 


dp. 


Here | 12 1 /2irt is our old friend from before and it remains to calculate asymp- 
totically (as t— >0) the integral J n e~ 28lt dp. To do this consider the curve T(8) of 
points in 12 whose “distance” from T is 5. (See Fig. 4.) 
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For small enough 5, T(5) is well defined (and even convex) and the major 
contribution to our integral comes from small 6. 

If L( 5) denotes the length of T(S) we have 


/ > /»5 0 

e~ 2h 2,1 dp — I e~~ 282/t L(5) dd + something less than | flj e“ 2, ^ 2/ * 
o. d n 


and hence, neglecting an exponentially small term (as well as terms of order t) 


f. 


e-^i'dp^Vt 


s. 


8 0 /^t 


/; 


L 


e 2x L(x\/t) dx ~ y/tL J e~ 2x dx — — ■ 

4 


where L — L{ 0) is the length of T. 
We are finally led to the formula 


Y e~ Xnt 


JYi 

2i rt 


L 1 
4 \/2'irt 


> for t — > 0, 


and so we can also a hear” the length of the circumference of the drum! 

The last asymptotic formula was proved only a few years ago by the Swedish 
mathematician Ake Pleijel [2] using an entirely different approach. 

It is worth remarking that we can now prove that if all the frequencies of a 
drum are equal to those of a circular drum then the drum must itself be circular. 
This follows at once from the classical isoperimetric inequality which states 
that L 2 ^47r| Q \ , with equality occurring only for a circle. 

By pitch alone one can thus determine whether a drum is circular or not! 


Kp) 



Fig. 5 


1L Can the heuristic argument again be made rigorous? Indeed it can. First, 
we use the inequality 
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which is simply a refinement of the one used previously, namely, P fi (p|p; t) 
^ \/2irt, and which can be proven the same way. 

Next we need a precise lower estimate for Pfip\ p ; t) and this is a little more 
difficult. We “inscribe” the rectangle P(p, t) as shown in Fig. 5, where hit), the 
height of the shaded segment, is to be determined a little later. 

Let the side of R along the base of the segment be b(t) and the other side 
be w{t). It should be clear from the picture that the y-axis bisects the sides of the 
rectangle which are parallel to the x-axis. 

Now consider Pr(p\ p; t). This notation is perhaps confusing since it suggests 
that we are dealing with a boundary value problem in which the boundary 
varies with time. This is not the case. What we have in mind is the following: 
fix t , find Pi 2 (i )(p \r; r) which is defined unambiguously, and finally set r = t. 
The result is P^(p|r; t). A convenient expression is 


Pr(p | p ; 0 


=— \i 

2ir t l Z 


» / r 2b 2 i r 2b u 1 yi\ 

?.r p L“T”’J- exp L“TV ,+ TjJj 


X 


2 w 2 


r 2w 2 / 8 VI 

n 2 

t 

— exp 

( n H 1 

L t \ w ) J 


where 8 = 8 — h(t) =\\q—p\\ —h(t). Now let h(t)=e\/t and, assuming that /(p) is 
actually tangent to the curve (which for a convex curve will happen with at most 
a denumerable number of exceptional points q ), we have 


and consequently 


.. b(t) b(t) 

lim = lim 

«-0 h(t) t — >0 eVt 




1 + o(l). 


This is not quite enough, however, and one needs the stronger estimate 




1 + (?(\/ 1) . 


This will surely be the case, for example, if the curvature exists at q, for this 
would imply that h(t)~b 2 (t) and the o(^/t) term above would then be an 
enormous overestimate. Very mild additional regularity conditions at nearly all 
points q would insure o{\/t). Without entering into a discussion of these condi- 
tions let us simply assume the boundary to be such as to guarantee at least 

o(Vt). 

Since w(t) remains bounded from below as t— >0, we also have 



= 1 — e~ u2lt + exponentially small terms. 
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We are now almost through. We write (cf. Fig. 6) 

e ~ Xnt = f -P»(p 1 p; 0 dp > f -Pr(p I p; 0 ^p 

h« 1 ^0(eV<) 

(i + (7(vo) r 


2 7T ^ 


J (1 — e 252/ * + exponentially small terms) dp. 

fi(eV t ) 



Fig. 6 


Except then for exponentially small terms and the factor 1 +p(V 0 in front 
we have the integral 

f (1 - dp 

J fi(eV t) 


which, as before, can be seen to be asymptotically 


£l(e\/t) | — \Z2irt, 


where one neglects terms of order t and exponentially small terms. Since asymp- 
totically 1 £l(e\/t) | S2 1 — Ley/t one can obtain the inequality 

g 1 (L + e') 1 

2irt 4 \/ 2irt 

where e' is related in a simple way to €. Since e' can be made arbitrarily small, 
the asymptotic formula 

a| l i 

27r^ 4 yj 2i r/ 

follows. 

12. If our overall strategy of attack on the problem is right we should be 
able to go on and for points very close to a smooth boundary replace the bound- 
ary locally by suitable circles of curvature. 


CO 





Y, e~ Xnt > 

W =1 
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A result of Pleijel suggests strongly that for a simply connected drum with 
a smooth boundary (i.e. without corners and with curvature existing at every 
point) one has 


” ft 

T. e Xni ^ 

n= i 2rrt 


L 1 1 

4 y/ hrt 6 


Unfortunately I am unable to obtain this, for the exasperating reason 
that I am unable to get a workable expression for P si ( / o|p; t) if ft is a circle. 

Rather than yield to despair over this sad state of affairs let me devote the 
remainder of the lecture to polygonal drums, i.e. drums whose boundaries are 
polygons. This study will show beyond the shadow of a doubt that the constant 
term in our asymptotic expansion owes its existence to the overall curvature of 
the boundary. 

13. Before I go on I need an expression for Ps«> 0 )(p|r; t) where S(d 0 ) is an 
infinite wedge of angle 6 0 . In other words Ps(o 0 ) is the solution of 

dP 1 

— = — V 2 P 
dt 2 


subject to the usual initial condition Ps«j„)(p|r; <)— >8(p — ?), /-> 0, and vanishing 
as r approaches a point on either side of the angle do. 

This is a very old, very classical, problem and if 9 0 = Tr/m, with m an integer, 
it can be solved by the familiar method of images. For m not an integer, Som- 
merfeld invented a method which, so to speak, extends the method of images 
to a Riemann surface. A little later, in 1899 to be precise, H. S. Carslaw gave 
a more elementary approach in which Ps(# o) (p|r; t) is represented by a suitable 
contour integral. Carslaw transforms the integral into an infinite series of Bessel 
functions but for our purposes it is best to resist the temptation of Bessel func- 
tions and to reduce the integral to a different form. I shall skip the details 
(though some are quite instructive) and simply reproduce the final result 

Set 


»(«) = (1/2t t) y exp 

6—a—'K<_2k6Q 
<0— a-}-7 r 


r 2 — 2rpcos(0 — a — 2 kdo) + p 2 ~ 
2 1 _ 



where the summation ^3 is extended over k’s satisfying the inequality under the 
summation sign and p = (p, a), f = (r, 6). 
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Then 


Ps(6 0 )(p | r; t) = v(a) - v(-a). 

Note that if d 0 = 7r/m, with m an integer, the complicated integral is out, since 
the factor in front of it, to wit sin tt 2 /6q = sin inn, is zero; what remains in the 
resulting expression for v(a) — v(— a) is a collection of terms easily identifiable 
with those obtained by the method of images. 

Let us now assume that t/ 2 <d 0 <ir and see what P#(p|p; t) is in this case. 
In the expression for v(a) when we set d = a the inequality under the sign 
becomes — ir <2kd 0 <ir and only k = 0 is allowed. In v(—a) the inequality is 
2a — 7r<2kdo<2a-\-T and what k's to take depends on a. 

We see that: 


0 < a < 0o — > only k = 0 is allowed, 


< a < 0 O , only k = 1 is allowed, 


but for 6 0 — 7t/2 <a <tt/2 both k = 0 and k = 1 are allowed. (See Fig. 7.) 

Let us now put r = p (so that p = r ) and write down in detail the expressions 
for Ps(6 0 )(p\p', t) in the three sectors. For 0 < a < do — 7t/2 


exp 


Ps(p | p; t) = 


(1 — cos 2a) 


2i rt 


^sin 


+ sin - 


exp 


2 71 ’t 

r 2~ 

r r 2 

— 

exp cosh y 

t _ 

r°° L t J 


^irdot 


7T 7 r 

cosh — y — cos - — 
do do 


dy 


For 7r / 2 <a <d 0 


exp 

r 2- 
t _ 

exp 

r 

r 2 

cosh y 


4 irdot 

J-oo (7 r 

a ) 

IT 2 



cosh < — 

y + 2 ici — > - 

- cos — 



(do 

9o) 

do 


r 

r 2 

1 



dy. 


Ps(p | p; 0 = 


1 


exp 


t 


(1 — cos 2(d 0 — a) 


2wt 2ir 

+ the same two integrals as above 



CAN ONE HEAR THE SHAPE OF A DRUM? 


21 


7T 

2 



and finally, for 6 0 — t/ 2 <a<ir/2 


Ps(p | p ; t) = 



r 2 1 

r 2 

exp 

(1 — cos 2a) exp 

— (1 — cos 2 (do — a)) 


2i rt 

2irt 


2tI 

+ again the same two integrals. 

We should recognize r 2 ( 1— cos 2a) (and f 2 (l — cos 2(6q — a:))) as being 2d 2 
where 8 is the distance from p to a side of the wedge. 



14. To simplify matters somewhat let me assume that the polygonal drum 
is convex and that every angle is obtuse. 

At each vertex we draw perpendiculars to the sides of the polygon thus ob- 
taining N shaded sectors (where N is the number of sides or vertices of our 
polygon). 
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Now let p be a point in 12. Stuff diffusing from p will either “see” the boundary 
as a straight line or, if p is near a vertex, as an infinite wedge. 

We may as well say that the boundary will appear to the diffusing particle 
as the nearest wedge, and that consequently we may replace P fi (p|p; t) by 
P s(8 0 ) (p | p; t ) , where S{S 0 ) is the wedge nearest to p. 

Now, each Ps(p\ P\ t) has \/2irt as a term and after integration over 12 this 
gives the principal term |l2|/27r2. Next, each P^(pjp; t) contains two compli- 
cated looking integrals which have to be integrated over the wedge. 

Fortunately, the second of these integrates out to 0, while the first yields, 
upon integration over S(9 0 ), 



dy 


(1 + cosh y) 




This is only the contribution of one wedge; to get the total contribution one 
must sum over all wedges. 

Thus the total contribution is 


--Z(sin-) f 

87T 6q V So/ J 


dy 

7 r 7 r 

(1 + cosh y)( cosh — y — cos — 
So 


Finally, if p is in the shaded sector of the wedge S(S 0 ) we get, on integrating 
over the sector, 


- r *,/•' 

d e 0 — tt/2 d n 



r 2 

exp 

(1 — cos 2a) 

_ t _ 

27 rt 


exp 


+ 


(1 — cos 2(00 — a)) 


2irt 


rdr 


2 — C °t ( So 
2 tt 


and the total contribution from the shaded sectors is —\\/2i r^o Q cot (0 O -tt/2). 
The remaining contribution is easily seen to be 




er 2 * i*d8 


L 

4 -\J 2t rt 


* +i iy C otb„-P). 

2t rt 2 tt Qq \ 2 / 


Finally, for a polygonal drum 
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00 


z 


e 


—\nt 


n= 1 


| 0| L 1 

2 t rt 4 sj 2irt 


— f 00 — , 

8 t r 7^ V d 0 J (1 + cosh 3;) (cos w/ do y — cos7r 2 /0 o ) 


with the understanding that each d 0 satisfies the inequality tt/2 <Oo<t. If the 
polygon has N sides, and if we let N — »<*> in such a way that each 0 O — then 
the constant term approaches 


+ 



dy 

(1 + cosh y ) 2 


1 

6 


This should strengthen our belief that for simply connected smooth drums the 
constant is universal and equal to J. 


15. What happens for multiply connected drums? 

If the drum as well as the holes are polygonal the answer is easily obtained. 
One only needs Ps(e 0 )(p|p; t) for d 0 satisfying the inequality t<0o<2it and this 
is easily gotten from the general formula quoted above. 

Near the holes the diffusing particles will “see” concave wedges but nothing 
will change in principle. 

If we let all polygons approach smooth curves it turns out the constant 
approaches (1 — r)£, where r is the number of holes. It is thus natural to conjec- 
ture that for a smooth drum with r smooth holes 


E 



L 

4 


■\Z2tt t 


+ (1 ~ r) i 


and that therefore one can “hear” the connectivity of the drum! 

One can, of course, speculate on whether in general one can hear the Euler- 
Poincare characteristic and raise all sorts of other interesting questions. 

As our study of the polygonal drum shows, the structure of the constant term 
is quite complex since it combines metric and topological features. Whether 
these can be properly disentangled remains to be seen. 

This is an expanded version of a lecture which was filmed under the auspices of the Committee 
on Educational Media of the Mathematical Association of America. 
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